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Non-invasive quantitative myocardial perfusion assessment 



01.07.2003 



This invention pertains in general to the field of image analysis and more 
particularly to image analysis of a series of non-invasively acquired cardiac images and even 
more particularly to the non-invasive assessment of perfusion of the myocardium by image 
analysis of a series of non-invasively acquired cardiac images. 
5 The myocardium (MC) is the heart's muscular wall. It contracts to pump blood 

out of the heart, then relaxes as the heart refills with returning blood. Blood perfusion of the 
myocardium ensures functioning of the myocardium. Reduced blood perfusion of the 
myocardium is a cardiovascular disease. The quantitative assessment of myocardial perfusion 
is thus of importance for the diagnosis of this cardiac disease. Insufficient perfusion of the 

10 myocardium may be the result of the (partial) occlusion of a coronary artery, which is a 

disease called ischemia. Cardiac Magnetic Resonance Imaging (MRI) is becoming one of the 
preferred techniques for the diagnosis of this kind of cardiac disease. Myocardial perfusion is 
today calculated from first-pass myocardial perfusion images acquired with ECG-triggered 
contrast-enhanced MRI. Myocardial time-intensity profiles are measured from this series of 

1 5 cardiac images. 

According to a method of quantification of myocardial perfusion, local 
myocardial time-intensity profiles are firstly measured from scans acquired when the heart 
was at rest and then when it was stressed. Perfusion is then quantified from these profiles. An 
example for a perfusion parameter thus derived is the so-called myocardial perfusion reserve 

20 index (MPRI), which has proven to be a good indicator of myocardial ischemia, as has been 
disclosed in Wilke at al., Radiology 1997; 204: pp. 373-384 and Al-Saadi et al., Circulation 
2000; 101: pp. 1379-1383. The myocardial perfusion reserve index (MPRI) is defined as the 
maximum upslope of the time-intensity profile at stress divided by that of the time-intensity 
profile at rest In order to normalize for the speed of injection of the contrast agent, the 

25 myocardial maximum up slopes are first divided by the maximum upslope of the time- 
intensity profile of the left ventricle (LV), when calculating myocardial perfusion parameters 
such as the MPRI. This time-intensity profile is usually called the LV input function. 

Clinical validation has shown that the above described method may fail due to 
the fact that often only an insufficient LV input function is available and that therefore the 
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maximum upslope of this LV function cannot be measured accurately. This occurs due to e.g. 
patient breathing or to non-linear behavior of the scanner at high image intensity levels, so 
that there is no linear relation between the increase in contrast agent in the LV and the 
increase in intensity in the perfusion images. As a result, the value of the calculated perfusion 
5 parameters, such as the MPPJ, may not be very reliable or can even be completely wrong. 

Therefore, a need exists for a reliable way of non-invasively quantitatively 
assessing myocardial perfusion independently of the relation between the increase in contrast 
agent in a ventricle and the increase in intensity in the myocardium in the perfusion images. 

10 01X6 object of the invention is to provide a reliable way of non-invasively 

quantitatively assessing myocardial perfusion. Another object is to provide a way of non- 
invasively quantitatively assessing myocardial perfusion independently of the relation 
between the increase in contrast agent in a ventricle and the increase in intensity in the 
myocardium in the perfusion images. 

15 l^e present invention overcomes the above-identified deficiencies in the art 

and solves at least the above identified problems by providing a computer-readable medium, 
a device and a method for the quantitative assessment of cardiac perfusion according to the 
appended patent claims. 

Cardiac perfusion is assessed by means of cardiac perfusion parameters, i.e. 

20 parameters that quantify the uptake of the contrast agent in the myocardium. 

The general solution according to the invention is based on the fact that any 
myocardium contains parts that are normally perfused because otherwise the heart muscle 
would not be able to perform any pump action at alL More particularly, the perfusion in these 
normally perfused parts is taken as a reference of the speed of contrast medium injection, 

25 instead of deriving this information from e.g. the LV time-intensity profile. 

According to aspects of the invention a computer-readable medium, a device 
and a method for the quantitative assessment of cardiac perfusion are disclosed. 

According to one aspect of the invention a computer-readable medium is 
provided, which stores a computer program for the non-invasive quantitative assessment of 

30 cardiac perfusion from a series of cardiac images, such as a 1 ^-pass perfusion image series. 
This computer program is in use processed by a computer and has a plurality of code 
segments, wherein a first code segment selects one or more segments with normal perfusion, 
such that cardiac perfusion parameters of the remaining image segments are based on a 
perfusion parameter of said image segment(s) with normal perfusion. According to 
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embodiments of the invention, further code segments divide the myocardium depicted on the 
image series into segments, a time-intensity profile for distribution of a contrast agent in the 
myocardium is determined for each of these segments, and a mavimnm upslope for each 
time-intensity profile of myocardium image segments is determined. Furthermore, a normal 
5 maximum upslope is derived and a relative maximum upslope is calculated for each segment. 

According to another aspect of the invention an apparatus is provided, wherein 
said apparatus comprises means for executing the above computer program. 

According to a further aspect of the invention, a method for the quantitative 
assessment of cardiac perfusion from a non-invasively captured series of cardiac images, 
10 such as a 1 ^-pass perfusion image series is provided. The method comprises the step of 

selecting at least one segment with normal perfusion, such that cardiac perfusion parameters 
of the remaining image segments are based on a perfusion parameter of said image 
segments) with normal perfusion. According to embodiments, the method further comprises 
dividing a myocardium depicted on the series of cardiac images into image segments. In a 
15 further step, a time-intensity profile for distribution of a contrast agent in said myocardium is 
determined for each of said image segments. Then a m^xim^rn upslope for each time- 
intensity profile is determined. Subsequently a normal maximum upslope is derived and a 
relative maximum upslope is calculated for each image segment 

The present invention has the advantage over the prior art that it reduces errors 
20 in cardiac perfusion measurement due to e.g. breathing artefacts and/or non-linear scanner 
behavior. 



Further objects, features and advantages of the invention will become apparent 
25 from the following description of embodiments of the present invention, reference being 
made to the accompanying drawings, in which: 

Fig. 1 is a perspective view schematically illustrating the term "short axis"; 

Fig. 2 is a graph illustrating the intensity progress as a function of time of 
different cardiac regions (the time-intensity profile); 
30 Fig. 3 A is an exemplary cardiac MR image showing cardiac segments with the 

relative maximum upslopes at rest; 

Fig. 3B is an exemplary cardiac MR image showing cardiac segments with the 
relative maximum upslopes at stress; 
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Fig. 3C is a schematic illustration of the thresholded MPRI corresponding to 
Figs. 3A and 3B; 

Fig. 4A and 4B are a schematic illustration of Fig. 3A and 3B respectively; 
Fig. 5 is a schematic diagram showing a computer-readable medium according 
to an embodiment of the invention; 

Fig. 6 is a schematic diagram showing a device according to an embodiment 
of the invention; and 

Fig. 7 is a flow chart illustrating a method for the assessment of cardiac 
perfusion according to an embodiment of the invention. 



10 



MR cardio perfusion scans are obtained by recording a number of short-axis 
slices during a period of 20 to 40 seconds. The definition of a short-axis slice is illustrated in 
Fig. 1, wherein a long axis 101 and a short-axis slice 102 are shown. Furthermore the 

15 schematic illustration of cardiac components 100 shows the left ventricle 105, the right 

ventricle 106, the myocardium 103, the epicardial contour 107 and the endocardial contour 
104. The acquisition of the cardio perfusion scans is ECG controlled so that each set of dices 
represents the same phase of the heart cycle, usually starting from end diastole. A few 
seconds after the beginning of the scan, a contrast agent is injected into the patient and 

20 subsequently the heart is imaged. The heart is also imaged in the period prior to the injection 
of the contrast agent, to acquire so called baseline images. These baseline images are 
sometimes used to correct for intensity inhomogeneity. The contrast agent shows brightly in 
the MR images. In a MR 1 "'-pass perfusion image series, it is observed that the contrast agent 
first enters the right ventricle (RV), then the left ventricle (LV) and finally the myocardium 

25 (MC). The intensity I as a result of the contrast agent as a function of time t of a segment in 
the RV, the LV and the MC is shown in graph 200 in Fig. 2. Firstly the contrast agent is 
spread in me RV, as shown by the intensity increase of the dotted curve 201, which 
represents the RV bolus. Then the contrast agent is spread in the LV, as shown by the 
intensity increase of the dashed curve 202, which represents the LV bolus. The spreading of 

30 the contrast agent in the MC, is shown by the solid curve 203 and the maximum upslope is 
indicated with the straight line 204. The myocardial perfusion reserve index (MPRI) is 
defined as the maximum upslope of the time-intensity profile at stress divided by that of the 
time-intensity profile at rest In order to normalize for the speed of injection of the contrast 
agent, the myocardial maximum upslopes are according to the prior art first divided by the 
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maximum upslope of the time-intensity profile of the left ventricle (LV). As mentioned 
above, this is not always reliable. An exemplary image illustrating this problem is shown in 
Figs. 3A and 3B. 

Since the signal in the myocardium may be very noisy, generally the 

5 myocardium is divided into segments and the intensity measurements are averaged over these 
segments. These segments generally coincide with areas of the myocardium, which are 
supplied with blood from a certain coronary artery. In this way, if a reduced perfusion is 
observed in a myocardial segment, it can be traced back to the supplying artery. Depending 
on the level of noise in the images, this myocardial segment consists of at least one image 

10 pixel. Cardiac images having a sufficiently low noise level, allow assessment of myocardial 
perfusion at the pixel level. 

Quantitative perfusion analysis is performed by using the normalized upslope 
of the myocardial signal enhancement to derive the myocardial perfusion parameters or 
indexes, such as the myocardial-perfusion reserve index. 

15 In an embodiment of the invention according to Fig. 5, a computer readable 

medium 500 stores a computer program 501 for the non-invasive quantitative assessment of 
cardiac perfusion from a cine series. This computer program is in use processed by a 
computer 50 and has a plurality of code segments 51-57, wherein a first code segment 51 
divides a myocardium depicted on the cine series into segments. A second code segment 52 

20 determines a time-intensity profile for distribution of a contrast agent in said myocardium for 
each of said segments. Then a third code segment 53 determines a maximum upslope for 
each time-intensity profile and a fourth code segment 54 selects one or more segments with 
normal perfusion. A fifth code segment 55 derives subsequently a normal mavim^m upslope 
and a sixth code segment 56 calculates a relative maximum upslope for each segment Finally 

25 a seventh code segment 57 calculates the myocardial perfusion reserve index (MPRI) for 
each segment from a ratio of relative maximum upslopes derived at rest and at stress. 

A further embodiment of the invention is illustrated in Fig. 6, wherein a 
workstation 600 is provided for executing the above computer program 501. 

In another embodiment of the invention according to Fig. 7, a method 700 for 

30 quantitative assessment of cardiac perfusion from cine series is shown. The method 700 
comprises the following steps and starts with the step 70 by dividing the myocardium into 
small segments. Then the myocardial time-intensity profile is measured in each segment 
according to the next step 71 . In the subsequent step 72 the maximum upslope is measured 
for each time-intensity profile of the segments. Subsequently one or more 'normal' segments 
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are automatically selected in step 73 and a 'normal' maximum upslope value is derived from 
the selected 'normal' segments in step 74. The maximum upslopes of all other segments are 
then expressed as a percentage of the normal maximum upslope, which henceforth will he 
called the 'relative maximum upslope', according to step 75. Finally, in step 76, the MPRI is 
5 calculated as the ratio of the relative maximum upslopes at stress and at rest. 

The normal segments), may he chosen according to criteria such as: 
the segment with the highest maximum upslope, 

the average of N segments with the N highest maximum upslopes, e.g. N = 3 
for a division of the myocardium into 24 segments, 
1 0 average of N segments with the N highest maximum upslopes, hut 

excluding those that have an upslope helow a certain threshold in order to avoid inclusion of 
ischemic segments. 

Using the first example, the highest upslope is taken as a reference, which 
means mat the relative upslopes of all other segments are always helow 100%. Furthermore, 

15 the MPRI is never much larger than 1.0. In fact, values around 1.0 indicate myocardial tissue 
that is well perfused hoth at rest and at stress, whereas significantly lower values indicate 
potentially ischemic tissue, which is not sufficiently perfused at stress. 

Figs 3A, 3B, 3C, 4A and 4B show an iUustrative example of performing an 
embodiment of the invention. Figs. 3A and 3B show two exemplary first-pass myocardial 

20 perfusion images, wherein the images are divided into image segments and the relative 
m axi mum upslopes at rest and at stress are respectively indicated in the images. Fig. 3A 
shows the maximum upslope at rest Fig. 3B shows the maximum upslope at stress. The 
magnitude of the maximum upslope is coded by means of different grey levels in me images, 
wherein black corresponds to low values of the maximum upslope and white corresponds to 

25 high values of maximum upslope. Greyscale values between black and white indicate 
maximum upslope values between these extreme values. Fig. 4A and 4B are schematic 
illustrations of Fig. 3A and 3B respectively. The cardiac images 300, 310, 320, 400, 410 are 
divided into exemplary twelve sectors 1-12 delineated by lines 401, 402. The myocardium is 
delineated by concentric contour lines 403, 404 and subdivided by a concentric line 405. In 

30 this way the myocardium depicted on the images is divided into segments 406, 407. For these 
exemplary segments, a perfusion parameter is calculated, and may be visualized by means of 
greyscale coding, as in Figs. 3A, 3B, 4A and 4B, or by means of other appropriate coding 
methods, such as color coding. 
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Fig. 3C shows a schematic illustration of the thresholded MPRI corresponding 
to Figs. 3A and 3B as well as Figs. 4A and 4B. In the example of Fig. 3C, a white segment 
indicates that the MPRI is above a certain, e.g. 75%, and a black segment indicates that the 
MPRI is below said threshold. 
5 These patient images cannot be analyzed in the conventional way due to the 

low quality of the LV input function. According to Fig. 3C, the MPRI calculated with the 
novel method clearly shows a large insufficiently perfused area at the bottom as well as two 
smaller areas on the right and left side. 

Applications and use of the above described computer readable medium, 
10 device, and method according to the invention are various and include exemplary fields such 
as detection of normal, ischemic and infarcted myocardial tissue using non-invasive medical 
imaging. 

The present invention has been described above with reference to specific 
embodiments. However, other embodiments than the preferred above are equally possible 

15 within the scope of the appended claims, e.g. different ways of determining 'normal' 

myocardial segments than those described above, performing the above method by hardware 
or software, etc. The invention is further not limited to a certain imaging technique. 
Generally the cine sequences used for analysis by the invention may be captured by any 
suitable medical imaging technique, i.e. any imaging technique for imaging inside parts of a 

20 mammal body providing suitable cardiac cine sequences. These imaging techniques comprise 
Magnetic Resonance (MR), Computer Tomography (CT), Nuclear Medicine (NM) and 
Ultrasound (US) Imaging. 

The described computer program, method and device are not limited to 
myocardial segments consisting of a certain number of pixels. Provided that the image 

25 quality is sufficient, the computer program, method and device work on pixel level. Instead 
of choosing "normal" segments, "normal" pixels, etc. are chosen. 

The present invention is also not limited to specific myocardial perfusion 
parameters, but may be used for derivation of any myocardial perfusion parameter or 
myocardial perfusion indexes. The present invention has been described above with reference 

30 to specific embodiments and calculation of a MPRI for myocardial assessment. An example 
of another relevant parameter, which may be assessed by the present invention, is the 
accumulated inflow. The accumulated inflow is sometimes also called the area under the 
time-intensity profile curve, which can be calculated by integrating the segmental time- 
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intensity profile in a certain time window, e.g. the period during which the contrast agent 
passes the myocardium for the first time. 

Furthermore, the term "comprises/comprising" when used in this specification 
does not exclude other elements or steps, the terms "a" and "an" do not exclude a plurality 
and a single processor or other units may fulfill the functions of several of the units or 
circuits recited in the claims. 



PHNL030839EPP 

9 01.07.2003 

CLAIMS; 



1 • A computer-readable medium (500) having embodied thereon a computer 

program (501) for the non-invasive quantitative assessment of cardiac perfusion from a series 
of cardiac images comprising image segments, said computer program being for processing 
by a computer (50), said computer program comprising 
5 a first code segment (73) selecting at least one image segment with normal 

perfusion, such that cardiac perfusion parameters of the remaining image segments are based 
on a perfusion parameter of said at least one image segment having normal perfusion. 

2. The computer-readable medium according to claim 1, comprising 

10 a second code segment (70) dividing a myocardium depicted on said cardiac 

image series into said image segments, 

a third code segment (71) deteimining a time-intensity profile for distribution 
of a contrast agent in said myocardium for each of said image segments, 

a fourth code segment (72) determining said perfusion parameter for each of 
15 said time-intensity profiles of said image segments, 

a fifth code segment (73) deriving a normal perfusion parameter from said at 
least one image segment having normal perfusion, and 

a sixth code segment (75) calculating a relative perfusion parameter for each 
of said segments with relation to said normal perfusion parameter. 

20 

3. The computer-readable medium according to claims 1 or 2, said computer 
program comprising a code segment (76) for calculating a ratio of cardiac perfusion 
parameters derived at stress and cardiac perfusion parameters derived at rest for each image 
segment. 

25 

4. The computer-readable medium according to claim 3, wherein said ratio of 
cardiac perfusion parameters is a myocardial perfusion reserve index (MPRI). 
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5 - ^ computer-readable medium according to claim 4, wherein said MPRI is 
calculated from relative maximum upslopes derived at rest and at stress. 

6 - The computer-readable medium according to claim 3, wherein said ratio of 

5 cardiac perfusion parameters is a thresholded MPRI being calculated by thresholding a ratio 
calculated from relative maximum upslopes derived at rest and at stress. 

7 - The computer-readable medium according to any of the preceding claims, 
wherein said perfusion parameter is used for visualizing insufficiently perfused myocardial 

10 areas comprising at least one of said image segments. 

8 - The computer-readable medium according to any of the preceding claims, said 
first code segment comprising selecting an image segment with the highest perfusion 
parameter value of all image segments as the image segment having normal perfusion, 

15 wherein a high perfusion parameter value is defined as good perfusion. 

9 - The computer-readable medium according to any of claims 1 to 7, said first 
code segment comprising selecting an average metric calculated from N image segments with 
the N highest perfusion parameter values , wherein N is an integer number significantly 

20 lower than the total number of image segments. 

10 - The computer-readable medium according to claim 9 comprising a code 
segment for thresholding said N highest perfusion parameter values. 

25 1 1 • The computer-readable medium according to any of the preceding claims, 

wherein the series of cardiac images is series of first-pass myocardial perfusion images. 

12 - The computer-readable medium according to any of the preceding claims, 

wherein said image segment is at least one image pixel. 



30 



13 - A device (600) being adapted for the quantitative assessment of cardiac 

perfusion, said apparatus comprising means for executing the computer program according to 



claim 1. 
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14. The device according to claim 13, said device being a workstation. 

15. A method for the quantitative assessment of cardiac perfusion from a non- 
invasively captured cardiac series of cardiac images comprising image segments, said method 

5 comprising selecting at least one image segment with a normal perfusion, such that cardiac 
perfusion parameters of the remaining image segments are based on a perfusion parameter of 
said image segment with normal perfusion. 

16. The method according to claim 15, comprising 

10 dividing a myocardium depicted on said cardiac image series into image 

segments, 

deteimining a time-intensity profile for distribution of a contrast agent in said 
myocardium for each of said image segments, 

determining said perfusion parameter for each time-intensity profile, 
1 5 deriving a normal perfusion parameter, and 

calculating a relative perfusion parameter with relation to said normal 
perfusion parameter for each of said image segments. 

17. The method according to claims 15 or 16, comprising calculating a ratio of 
20 myocardial perfusion parameters derived at stress and myocardial perfusion parameters 

derived at rest for each image segment 

18. The method according to claim 17 wherein said ratio is a MPRI calculated as a 
ratio from relative maximum upslopes derived at rest and at stress. 

25 

19. The method according to claims 16 or 17, wherein said perfusion parameter is 
a maximum upslope of a time-intensity for distribution of a contrast agent in a myocardium 
comprising said image segments. 

30 20. The method according to claims 1 5 to 19, wherein said assessment is used for 

the diagnosis of cardiac diseases. 



21 . The use of a workstation for performing the method according to any of 

claims 15 to 20. 
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ABSTRACT: 



A computer-readable medium, a device and a method for the quantitative 
assessment of cardiac perfusion. A myocardium depicted on a series of cardiac images is 
divided into image segments, wherein the image segment comprises at least one image pixel- 
Then a cardiac perfusion parameter is determined for each of said image segments. Then at 
least one image segment with a normal perfusion parameter value is selected. Subsequently 
cardiac perfusion parameters of the remaining image segments are based on said normal 
perfusion parameter value of said image segment with normal perfusion. According to an 
embodiment, the perfusion parameter is a maximum upslope of a time-intensity profile for 
distribution of a contrast agent in said myocardium. A normal maximum upslope is derived 
for at least one image segment and a relative maximum upslope is calculated for each 
segment with relation to the normal maximum upslope. Based on these values, a ratio of 
myocardial perfusion parameters derived at stress and myocardial perfusion parameters 
derived at rest for each segment may be calculated. For example a myocardial perfusion 
reserve index (MPRI) is calculated for each image segment as a ratio of the relative 
maximum upslopes derived at rest and at stress. 
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